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Abstract. We discuss the so-called two-temperature model in linear thermoe¬ 
lasticity and provide a Hilbert space framework for proving well-posedness of the 
equations under consideration. With the abstract perspective of evolutionary equa¬ 
tions, the two-temperature model turns out to be a coupled system of the elastic 
equations and an abstract ode. Following this line of reasoning, we propose an¬ 
other model being entirely an abstract ode. We highlight also an alternative way 
for a two-temperatnre model, which might be of independent interest. 
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0 Introduction 

Chen and Curtin jT] and Chen et al. mm have given the formulation of the theory 
of heat conduction related to a deformable body which is based on two different 
temperatures. Here the first one is the conductive temperature, (j) and the other 
one is the thermodynamic temperature, 9. Chen et al. [2] discussed that these two 
temperatures are equal in absence of heat supply in the case of time independent 
situations and the difference between these two temperatures is proportional to 
the heat supply, where, as generally in the cases of time dependency, these two 
temperatures are different. Before these studies, by doing the study of transient 
coupled thermoelastic boundary value problem in half space, Boley and Tolins 
H] gave the conclusion that the strain and two temperatures are found to have 
explanation in the form of a wave plus a response taking place immediately through 
the body. The uniqueness and reciprocity theorems for the two- temperature 
thermoelasticity theory in case of a homogeneous and isotropic solid was reported 
by Lesan [5]. Subsequently, investigations were carried out on the basis of this 
theory by several researchers like, Warren and Chen Warren [7], Amos [8] etc. 

This theory (2TT) has drawn attention of researchers in recent years and some 
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specific features of this theory are reported (see jH (TOl HH HSJ [T3l [TH [151 ESI, 
and the references there-in). 

A structural formulation for linear material laws in classical mathematical physics 
was introduced by Picard |20) who considered a class of evolutionary problems 
which covers a number of initial boundary value problems of classical mathemat¬ 
ical physics. The corresponding solution theory is also established in |20) . Prior 
to this, Picard |21| also reported the structural formulation for linear thermo¬ 
elasticity in nonsmooth media. Recently, Mukhopadhyay et al. |22) have studied 
various models of thermoelasticity theory and have shown that the these models 
can be treated within the common structural framework of evolutionary equations 
and considering the flexibility of the structural perspective they obtained well- 
posedness results for a large class of generalized models allowing for more general 
material properties such as anisotropies, inhomogeneities, etc. It should be noted 
that evolutionary equations in the form just discussed have also been studied with 
regards to homogenization theory, see e.g. [251 [2n 123) . The aim of this note is 
to analyze the two-temperature thermoelastic model given by Chen and Curtin 
[T] as a first order system within the framework of evolutionary equations, see 
e.g. |26| . The models of thermoelasticity we shall discuss have been originally con¬ 
ceived as constant coefficient models. There is little harm in this assumption at 
this point, since we shall dispose of this simplihcation completely, when we discuss 
more general models in the last section. An alternative two-temperature thermoe¬ 
lastic model is proposed in which we can avoid involving roots of an unbounded 
operator. It is believed that the general perspective on two-temperature thermoe¬ 
lasticity to be presented may shed some new light on the theory of homogenization 
of such models. 

In sectionlH we discuss the functional analytic background needed, for discussing 
the two-temperature model. Section [2] discusses the two-temperature model in de¬ 
tail. In this section, we will also give a suitable Hilbert space framework allowing 
for well-posedness of the respective equation. An observation in section [2] is that 
the heat equation part is replaced by an abstract ode with infinite dimensional state 
space. More precisely, in the heat equation part the only unbounded operator in¬ 
volved is the time-derivative. Having realized this property of the two-temperature 
model, we propose in the two concluding sections [3] and H) alternative systems of 
thermoelasticity. The first one being entirely an abstract ode in the sense just 
discussed. The second one describes a possible alternative model, which does not 
involve square roots of operators. 


1 Functional Analytic Preliminaries 

In this section, we shall elaborate on some standard concepts in functional analysis 
needed in the following. Most frequently, we will have occasion to use the square 
root and the modulus of an operator: 

Definition 1. Let Hq, Hi be Hilbert spaces. Let C: D{C) C LTq Hq be 
a non-negative definite, selfadjoint operator, that is, for all (p E D{C) we have 
{(j),C(j)) > 0 and C = C*, then y/C is defined as the unique non-negative definite 
operator satisfying yfCyfC = C. For A\ D(A) C Ho Hi, a closed and densely 
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defined linear operator, we define the modulus of A, |A|, by 


1^1 := VA*A. 

Recall that D{A) = DdAI) and that ||A(/>|| = || |^| (f>\\ for all (j) G D{A). We 
record the following standard fact: 

Proposition 2. Let Hq, Hi be Hilbert spaees, A: D{A) C Hq —>• Hi densely 
defined, closed, linear. Then 


with 


A = a(^iJi + \A\^^ eL{Ho,Hi). 

aUi + \A[ 


< 1 . 


Proof. Let A = U \ A\ with a partial isometry U, being in particular a contraction, 


-1 


we have by the spectral theorem, 1 + G D(|A|) for all (j) G Hq, and 

thus, 

0 = |A| 

<\\(j)\\, {(j)GHo), 

establishing the boundedness and the norm-estimate of the operator A 1 -|- |^| 

As the operator 1 -|- A is densely defined, for the asserted equality in 

the proposition, it suffices to establish the inclusion 


-1 


\Jl + \A* 


-1 


AC A 




-1 


( 1 . 1 ) 


Next, we prove (HH): For this, by induction, we show the inclusion 

{1 + AA*)-'^ACA{1 + A*A)-'^ (nGN). (1.2) 

For proving the latter inclusion for n = 1, observe that for G D{AA*A), we have 

(1 + AA*)Afi = A{1 + A*A)cl). 

Hence, substituting '.= [1 + A*A)(f), we get 

A{1 + A*A)-^fi = (1 + AA*)-^Afi. 


So, for every nGN the inductive step can be shown as follows: 

(1 + = (1 + AA*)-^(1 + AA*)-^A 

C {1 + AA*)-'^A{1 + A*A)-^ 
C A{1 + A*A)-^{1 + A*A)-^ 
= A{1 + A* A)-^^+^\ 
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For the proof of m, we recall that for every real number x > 0, with |x| < 1 the 
binomial series gives 

n=0 ^ ' 

Putting Xg := —ey{l + y)~^ for some y > 0 and e g]0, 1[, we have \xe\ < e,l + Xe = 
(1 + (1 — e)y) (1 + which also leads to 

y/l + xe = \J+ £)y) (1 + y)~^ -t \l{l + y)~^. (1.4) 

Moreover, plugging Xg into the series , we arrive at 

E = E (T) 

n=0 ^ ' n=0 ^ ^ 

= E(‘f) {(-^rs-d+ !/)-”). 

n=0 ^ ^ 

By the functional calculus for selfadjoint operators, we may replace in the latter 
expression y by A*A and AA*, respectively. Thus, for e g]0, 1[, we set 

Bi,e ■■= f; ) ((-£)" {AA*) - (1 + {AA*))-^) , 

= E {A-Ar{l + (A-A ))-’'). (L5) 

Note that and i? 2 ,e define bounded linear operators. Moreover, by the spectral 
theorem (write AA* and A* A as multiplication operators in a suitable L^-space), 
we get invoking (|1.4I) 

^J{l + AA*)~^ and B 2 ,e ^ \l{l + A*A)~^ (1.6) 

as e —)• 1 in the strong operator topology. Thus, for e g]0, 1[ we get with the help 
of (II.2p and (II.5^ : 

Bi^^A = ((-£)" (AA*) ^ (1 + (AA*))-^) A 

n=0 '' ' 

c ^ {{-ef {AA*) ^A (1 + {A*A))-^) 

n=0 k / 

= ((-e)" A (A*A)^ (1 + (A*A))-^) 

= ^fl (^M)” (1 + {A*A))-^) 

n=0 V ^ / 

= AB2,e. 

Thus, the closedness of A, together with (II.6p yields the asserted inclusion (11.11) . 

□ 
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Another fact used in the following is mentioned in the next proposition. 
Proposition 3. Let Hq, Hi Hilbert spaces, A: D{A) C Hq — )• Hi densely defined, 
closed, linear, k G L{Hi) with 0 G Then kA is densely defined and closed 

and we have 

(nA)* =A*n*. 

Proof. The operator kA is clearly densely defined. Moreover, if {4>n)n is a sequence 
in D{A) such that {4>n)n and {KAfin)n are convergent to ip G Hq and r] G Hi, we 
infer, by the continuous invertibility of k and the closedness ot A, ip ^ D{A) and 
Afi = K~^'ip. Hence, kA is closed. The equality (kA)* = A*k* is also easy. □ 

Next, we briefly recall the functional analytic setting in which we are going to 
discuss the two-temperature model later on. A more detailed discussion can be 
found in |20l [26] or (particularly concerning the time-derivative) in [27] . See also 

m 

Definition 4. Let n > 0, H Hilbert space. Define (M, H) to be the space 
of (equivalence classes of) square integrable functions /: M —)■ 77 with respect 
to the measure with Lebesgue density x i—>■ Denote the space of L^- 

functions / with distributional derivative f representable as (M, i7)-function 
by Hi,^i (M, 77). Define 

do : 77,,1 (M, 77) C (M, 77) ^ (M, 77), / ^ /'. 

Note that we will not notationally distinguish between the time-derivative re¬ 
alized as an operator in L^(M, 77i) and L^(M, 772) for possibly different Hilbert 
spaces 77i and 772. The reason of introducing this particularly weighted L^-space 
is the fact that do becomes a continuously invertible operator. In fact, one has 
||5(("^|| < 1/j^, see [27] . 

For a closed and densely defined linear operator C: D{C) C Hq —>• 77i between 
the Hilbert spaces Hq and 77i, the lifted operator as an abstract multiplication 
operator from (M, 77o) to L^(M, 77i) will be denoted by the same notation. 
With these conventions, we can come to (a special case of) the solution theory first 
established in [20] . We mention here possible generalizations to non-autonomous 
mm) or non-linear frameworks (laiisQi). Denoting the range of an operator 
Mq by R{Mq) and its kernel by N {Mq) we recall the following general solution 
theory result from |20[ |26] . 

Theorem 5. Let 77 Hilbert space, Mq = Mq,Mi G L{H), A: D{A) C H ^ H 
skew-self adjoint. Assume there exists c > 0 such that {Mo<p, (p) > c{(p, (p) and 
‘Tit {Mi'ip,^) > c{ip,'ip) for all (p G R{Mo), ip G N{Mo). Then there exists vq P 0 
such that for all v > vq the operator sum 


B := doMo + Mi+A 

is closable as an operator in Lip (R, 77) and the closure B is continuously invertible 
in Lip (R, 77). Moreover, B is causal in the sense that given f £ L^ (R, 77) with 
the property that f = 0 on (—oo, a] for some a G R, then B f = 0 on (—oo, a]. 

The latter theorem tells us that the non-homogeneous problem Bu = / admits 
a solution for all f £ L‘p (R, 77) given v sufficiently large. In [26] it has been 
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shown how to invoke initial value problems in this context. Note that it is also 
possible to show that the solution u does not depend on the parameter v, that 
is, let > 0 be sufficiently large then the solution operators and 

established in (M, i?) and i?), respectively, coincide on the intersection 

of the respective domain, that is, on (M, H) n (M, H). 

Later on, we will also need the following operations skew : —>■ 

I (a — and sym : ^4 i-)- ^ (A + . 


2 The two-temperature model 


In this section, we shall have a deeper look into the two-temperature model found 
in [1]. For this, however, we have to introduce several vector analytical operators. 
In the whole section, we assume we are given an open set II C M"'. 

O 

Definition 6. We denote by Coo {A) the set of smooth functions with compact 
support. Then, we define, as usual. Grad to be the symmetric part of the 3x3- 
matrix-valued derivative of a smooth vector field <h, grad(/) be the gradient of a 
smooth function (f) and Div T and div ^ be the row-wise and the usual divergence 
for a smooth matrix-valued functions T and a smooth vector-valued function ip, re¬ 
spectively. Reusing the notation Grad, grad, Div and div for the respective L^(D)- 


o 

realizations, we further define Grad := Grad 

O 

o , Div := Div 



Coo(£^)^ sym 



grad := grad o , div := div 

Ooo 

joints — Div := ( Grad 


Coc(0)'= 


grad 


Coo(n) 


V Coc(n)' 

, —grad := (div 


Cooiny 


and their respective (Il)-type ad- 

, — div : = 

Note that here Div maps from 


1 Div 


' O o o' 

V 

sym 



and Grad maps into the Hilbert space (D) := (D,sym of 3 x 3- 

symmetric-matrix valued L^-type mappings. 

In the so-called two-temperature models of Chen and Gurtin |T] , apart from the 
temperature 9 another temperature (p, the conductive temperature, is introduced 
(together with a reference temperature Tq G ]0, oo[) such that 


9 — {(p — To) = —a div q. 


( 2 . 1 ) 


Here a G ]0,oo[ is a parameter, called the two-temperature parameter. Assuming 
homogeneous Dirichlet boundary conditions, Fourier’s law is then formulated in 
terms of the conductive temperature as 

q = -Kgrad {(p - Tq) , (2.2) 

where k G is a selfadjoint operator with k > c > 0. In addition, 

the two-temperature system consists of the heat equation with mass density ^>o £ 
L°°(D), Qq'> cq> 0, that is. 


do {goTorj) -F div g = qqQ 

















or - for our purposes - more conveniently 

do (gori) + div (q/To) = QoQ/To, (2.3) 

where q is the heat flux as in (12.2p . rj is the entropy and Q is the heat source. For 
the entropy rj we have the material law relating the entropy to the temperature 9 
and the strain tensor £ = Gradn, u the displacement, 

eoTov = qoX9 + To'j*£ (2.4) 

for some scalar A > 0, an operator 7 G L{L‘^{Q), Next, the strain tensor 

£ = Gradu is related to the stress tensor a and the temperature via the elasticity 
tensor C = C* G L (Lgyjj^(Q)) being strictly positive definite and 7 in the following 
way 

£ = C-^a + C-\9. (2.5) 

The two-temperature model is completed by the balance of momentum 

— Div cr = (2-6) 

for some given external force F. 

In the following, we will show that Theorem [5] is applicable to the equations 
dsn, dOI), Q, (EH, (EH, EH- Hence, the Hilbert space setting introduced 
in the previous section provides a functional analytic framework such that for all 
right hand sides F and Q there exists a unique solution to the two-temperature 
model depending continuously on F and Q. So, the task to be solved in the next 
lines is to find the right unknowns and, hence, the right operators Mq, Mi and A 
making Theorem [5] applicable. 

It should be noted that our reformulation of the two-temperature model reveals 
that the introduction of the second temperature transforms the heat equation into 
an ordinary differential equation with infinite dimensional state space. 

A first step towards our main goal in this section is the following observation 
yielded by ( 12 .ip and ( 12 . 211 : 

Proposition 7. Let k = k* G be strictly positive definite. Assume that 

Tg ,a G ]0,oo[ and 9 G and q G D(div), (j) G iA(grad) satisfy (12.111 and 

(|2.2p . Then TOf/Q Ka ■= ^fotK^fa we have 



Proof. Plugging in Fourier’s law we can rewrite (12.ip as 

0=^1 — div Kograd^ — Tq) . (2.8) 

The operator y/Iq^grad : iA(grad) C L?‘{£L) —)• is a closed densely defined 

linear operator, since k and hence, y/Ffi are boundedly invertible, see Proposi¬ 
tion [3j Moreover, its adjoint is given by ^grad^ y/Ffi = — diVy^AT^ (Proposition 

^Of course here yfoLKyfa = an, but we prefer to write it in this more symmetric fashion, since in the eventual 
first order model equations a can be chosen more generally, i.e. as a continuous, selfadjoint, strictly positive 
definite operator, without affecting well-posedness. Also k will be allowed to be a continuous, selfadjoint, 
strictly positive definite operator. 
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[3]) and thus, — divK^grad is a selfadjoint, non-negative operator. In particular, 
1 — divK^grad is boundedly invertible. Hence, rephrasing (12.2D in terms of the 
temperature 6 we are led to 


q = —Kgrad — div Kagradj — div K^gradj (</> — Tq) 

= — Kgrad — div Kq, grad') 6. 


Next, applying Proposition [2] to A := y/^grad we infer 



and 



yielding the assertion. 


□ 


With the latter observation, we are in the position of rewriting the two-temperature 
model as a system in the spirit of Theorem (5) 

Theorem 8 . Let k = k* E C = C* e , 7 E L(L2y^(Sd), 

^>0 € \,a,To E ]0,oo[. Then the system (12.ID . (12.2p . (|2.3D . (I 2 '4p . (|2.5D . 

(1^ may be rewritten into 


{doMo + Mi+A)U = J 


(2.9) 


with Oqu = V and 

( 

U = 


V 

a 


\ 


\ 1 - v^grad div ^q/Tq J 

where ^/an^/a = Ka, 


Mo = 


o\ 

0 


/ £>0 0 0 

0 C-i C-I7 

0 7*C-1 + 0 

V 0 0 0 0/ 


j = 


( QoF \ 
0 

qoQ/Tq 

V 0 


A = 


f 0 - Div 0 0 

—Grad 0 0 0 

0 0 0 0 

Vo 0 0 0 / 
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and 


Ml 


^1,32 


/O 0 

0 


0 0 

0 

0 

0 0 

0 

“-^1,32 


\0 0 Mi_32 To / 



= v^grady 1 — div K^grad . 

In particular, there exists I'o > 0 such that for all i' > vq the equation in (1^ 
admits for every J ^ (M, 0 © L?‘{Q) 0 L^(0)^) a unique so¬ 
lution U e D {doMo + Mi + A) C Ll (M, © L‘^{n) © L‘^{Qf). 


The solution operator is continuous and causal. 


Proof Before computing that the equation (doMo + Mi A)U = J is a reformu¬ 
lation of the two-temperature model, we establish the well-posedness issue first. 
For this, note that Mq = Mg and A = —A*. Next, we check that Mq is strictly 
positive definite on its range. For the purpose of symmetric Gauss elimination, we 
define the transformation matrix 


/I 0 0 0\ 
0 1 7 0 
0 0 10 
\0 0 0 1 / 


Hence, 



/I 0 0 0\ 


/I 0 0 0\ 

1 _ 

0 1 -7 0 

s* — 

0 10 0 


0 0 10 

? ^ — 

0 7* 1 0 


Vo 0 0 i) 


Vo 0 0 1/ 


5-1 = 


We compute that 

(5-i)*Mo5-i 


/I 

0 

0 

0 

/ 

0 

1 

0 0 


0 

- 7 * 

1 0 


VO 

0 

0 1) 

V 


Qo 0 

0 G-i 
0 7*C-i 
0 0 


(5-1)* = 


0 

c-^1 

)oio ^A + 7*C 

n 




1-1 


0 

-I7) 0 


fl 0 0 0 \ 
0 10 0 
0 -7* 1 0 

Vo 0 01 / 


\ 

/lO 

0 

0\ 


0 1 

-7 

0 


0 0 

1 

0 

/ 

0 

0 

0 

1/ 


/ go 0 0 0\ 

0 c-i 0 0 

0 0 Qo'^o 0 

\ 0 0 0 0 / 


Next, as bijective transformation 5 reduces the space 

R := L^{PLf © © l2(H) © {0}, 
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we infer R{Mq) = R. Moreover, for </> G ii we compute 


= {MoS-^S<l>,S-^S4>) 

= {{S-y MoS-^Sct>,S<t>) 

> c((/>, 0) 


for some c > 0. On N(Mo), the operator iHeMi, the real part of Mi, is given by 
multiplication by Tq > 0. Hence, the assertion concerning well-posedness follows, 
once we have established that Mi defines a bounded linear operator. This, however, 
is a direct consequence of Proposition [2l Indeed, 


y/Kgia,A\J 1 — div Kograd (j) 


v^grad 


—1 


1 + 


v^grad 




As a next step we proceed in showing that the two-temperature model admits 
the asserted reformulation. For this, in turn, it suffices to observe the following 
consequence of the equations (12.4p and (12.511 : 


QoToT] = QoXe + To-f*S 

= goXe + Toy* (C"V + 0-^0) . 


Hence, 

Moreover, from £ 


0oV = {QoTQ^X + -f*C + V. 
Gradu and d^u = u it follows that 

doS — Gradu = 0 


and the balance of momentum (j2.6l) reads as 


godov — Div a = goF. 

Recalling (12.7p from Proposition [71 we note that 


div (ij/To) = 





which eventually establishes the assertion. 

□ 

Note that Mi ^32 has moved from its place in A for the limit case a = 0 to the 
material law. 
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Remark 9. Symbolizing non-vanishing entries in the block operator matrices under 
consideration by clearly, the pattern of Mq is 


Mo 


/★ 0 0 0 \ 

0 ★ ★ 0 
0 ★ ★ 0 

\ 0 0 0 0 / 


But the pattern of Mi is 


IHeMi = 


/O 0 0 0 \ 
0 0 0 0 
0 0 0 0 
Vo 0 0 ★/ 


, 3m Ml 


/O 0 0 0 \ 
0 0 0 0 
ooo/k- 
Vo 0 ★ 0 / 


Moreover, 

/ 0 - Div 0 0 \ 

. _ —Grad 0 0 0 

0 0 0 0 

Vo 0 0 0/ 


We see that the system has partly been turned into an ode in an infinite dimensional 
state space. 


3 A Two-Temperature, Two-Strain Model. 

In this section, we shall elaborate briefly on the possibility of developing an alterna¬ 
tive model, such that the whole pde-part in the two-temperature model discussed 
in the previous section vanishes. We start with basically the same model as in 
Theorem [HI As a preparation for deriving the two-temperature, two-strain model, 
we consider first the following system, which is unitarily congruent to the one in 
Theorem (8) 



/ ^'0 

0 

0 


( 

0 

1 


0 


0 



0 


VO 

0 

0 

0/ 


/ 



V 

e 


y^^grad div 


1/2 


^q/To 


+ 


/O 0 0 

0 

\ 

0 0 0 

0 


0 0 0 

“-^1,32 


yO 0 Mi^32 

To 

/ 


/ 

V 


V 

9 


\ 


1 — V^grad div 



^q/To 


+ 
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+ A 


/ 

V 


V 

C-^/^a 

e 


\ 


1 — \/^grad div 



1/2 


Vk ^q/To j 


( ^ \ 

0 

qoQ/Tq 

\ 0 


with 

/ 0 -DivC^/^ooX 

^ ^ -C^/^Grad 0 0 0 

0 0 0 0 ’ 

\ 0 0 0 0 / 

where 

Ai^i ,32 = v^grady^ 1 — div K^grad . 

Taking this as a starting point and substituting Cj^ := y/^Cy/^ for some /3 > 0, 
we may propose analogously a similar modification of the elastic part yielding 



( Qo 

0 

0 0\ 

( 

0 

1 

0 

1 

0 


0 

1 

* 

^To-'A + y^C-iy) 0 


^0 

0 

0 0/ 

( 


V 


^1 — -yC^Grad Div y/Cp yfC ^ a 


y \J\ — y/ vTiK-yagrad div y/y/an 
( 0 -M/io 0 0 0 \ / 


ay/K ^q/To j 


Ml 


,10 


0 0 0 


0 0 0 0 0 

0 0 0 0 -M /32 


V 0 

( ^ \ 

0 

qoQ/Tq 

\ 0 ; 


Grad Div y/C^ yfC ^ a 


0 h To / \ yl - ^/To / 


where now 


and 


Mi^io = —VCgv&dyJ 1 — Y^^Grad Div 


,^ 1,32 = v^grady /1 — div K^grad 


—1 
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Clearly, the pattern of Mq is still 


Mo 


/★ 0 0 0 \ 

0 ★ ★ 0 
0 ★ ★ 0 

Vo 0 0 0/ 


But the pattern of Mi is now 


symMi 


/O 0 0 0 \ 
0 0 0 0 
0 0 0 0 
Vo 0 0 ★/ 


skewMi 


/ 0 ★ 0 0 \ 
★ 0 0 0 
0 0 0 ★ 
Vo 0 ★ 0 / 


Moreover, 

/O 0 0 0\ 
0 0 0 0 
0 0 0 0 
VO 0 0 0/ 


We see that the system has completely been turned into an abstract ode. 
Remark 10. 


• Taking the general perspective used here into account, more complex ma¬ 
terials, for instance, the Maxwell-Cattaneo-Vernotte (MCV) model of heat 
conduction IMlISlEg, are easily possible to include all these models as in¬ 
troduced in HU. For implementing the MCV model we merely have to take 
Mq with the pattern 


Mo 


/★ 0 0 0 \ 

0 ★ ★ 0 
0 ★ ★ 0 

Vo 0 0★/ 


as strictly positive definite. 

• Moreover, if we change the parameter a (and /?) to be bounded, selfadjoint, 
strictly positive definite operators in an appropriate Hilbert spaces, we gain 
further flexibility for material modeling within the framework of the first order 
system. 

• Given the intricate rationale used in deriving the model in the first place 
it is somewhat disappointing to see that it merely serves to approximate a 
pde by an ode, which of course is always possible, compare e.g. the Yosida 
approximation or the above strategy, which amounts to replacing an un¬ 
bounded skew-selfadjoint operator A by the bounded skew-selfadjoint opera¬ 
tor Ay/l — aA^ ^ = y/l — aA^ ^A, a G ]0, oo[. 


4 An alternative two-temperature model 

In this section, we will make an attempt to establish an alternative two-temperature 
model from a purely structural point of view. For this, we proceed as follows: 
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Note that a transition to the ode setting can also be achieved for example by 


approximating A with A{1 + eA) (Yosida approximation). Indeed 

A{l + eA)-^^^^ A 


(4.1) 


point-wise on D (^). 

This way the occurrence of a square root (of inverses) of unbounded operators 
can be avoided. We assume the conditions of Theorem [H] For notational con¬ 
venience we set D := y^grad. Applying the idea of using (14.ip to our initial 
two-temperature model yields 


5, 


0 \ 

0 


/ QO 0 0 

0 C-i C-H 

0 (poTo-'A + y^C-^) 0 

\0 0 0 0 / 

/ u \ 

"" 0 ) + 

V (1 + e^DD*) ^~\lTo + eD0 / 

/O 0 0 0 \ 

/ 0 0 0 0 

0 0 eD*D{l+e‘^D*D)~^ -D* {l + DD*)~^ 

^0 0 D{l + s‘^D*D)~^ eDD* {l + e‘^DD*)~^+ToJ 

/ V \ 

I). 

-1 


V (1 T e'^DD*) y/E~qlTo + eD9 

( V 


-h A 




a 


U 


V (1 + e^DD*) ^~\/To + eDO 

( QoF \ 

0 

qoQ/Tq 

\ 0 / 


with 


A = 


( 0 - Div 0 0 \ 

—Grad 0 0 0 

0 0 0 0 

Vo 0 0 0/ 


This reduces indeed to 


doQQV - Diver = / 
do (cr -|- 70 ) — CGradu = 0 

do + (poTo-'A + 7*C^"S) O) - = QoQ/To 
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and finally 


D{l+e^D*D)-^d+ 

+eDD* (1 + e^DD*)~^ (^(1 + s‘^DD*)^~\/To + eDO^ + 

+ro ((1 + e^DD*)^~\/To + sDe'^ = 0, 

{eDD* + To + Toe^DD*) ^~^q/TQ+ 

+ (^(1 + e^DD*)-^ + £‘^DD* (l + e^DD*)~^ + eTo) D9 = 0, 

(1 + £“^00*) + £DD*^~\/To + (1 + £To) do = 0 , 

(1 + £^DD*) VK~^q + D (^£D*^~^q/To + 9 + = 0. 

Thus, 

^/^~\+{1 + £‘^DD*)~^ D(£D*^/^~\/To +9+ £To9^ =0, 

+ T» ((1 + e‘^D*D)~^ £D*^~\/To + (1 + eTo) (l + 0 ) = 0, 

(4.2) 


and hence, defining 

:= (1 + eTo) {l + £‘^D*Dy^9 + £{l + £^D*Dy^ D*^~\/To + To 

and recalling that D = -y/^grad, we end up with 

0 = (1 + £To)-^ (1 + £^D*D) (cj) -To)-^{l + eTo)-^ D*y^-y 

-to 

= (1 + eTo)“^ (1 - e^divKgrad) {(j) - Tq) + ^ (1 + eTo)“Mivg 
and (14.2p gives the Fourier law 


q + Kgrad 0 = 0 . 

Thus, using —divKgrad0 = divq, we get that 

0 = (1 TeTo)"^ (</>- To) + (1 + eTo)"^e^ divg + (1 + eTo)"^;^ divg 

-to 

= (1 + eTo)“^ {(j) -To) + ^ ((1 + eTo)“^eTodivg + (1 + eTo)“^ divg) 

-to 

= 4> — Tq + (div q - 

^ " To V 1 + £ToV 

This can also be written as 


0 - 1 - 


eTo 


1 + eTo 


((/> - To) = — div O' 


(4.3) 


Equation (14.311 represents the hnal relation satisfied by the two temperatures. 
The parameter e would be an alternative two-temperature parameter. 


17 




References 


[1] P. J. Chen & M. E. Gurtin , On a theory of heat conduction involving two 
temperatures, Z. Angew. Math. Phys. 1968;19: 614-627. 

[2] P. J. Chen, M. E. Gurtin & W. O. Williams, A note on non-simple heat 
conduction, Z. Angew. Math. Phys. 1968; 19: 969-970. 

[3] P. J. Chen, M. E. Gurtin & W. O. Williams, On the thermodynamics of non¬ 
simple elastic materials with two temperatures, Z. Angew. Math. Phys. 1969; 
20: 107-112. 

[4] B.A. Boley, I.S. Tolins, Transient coupled thermoelastic boundary value prob¬ 
lems in the half space, J. Appl. Mech. 1962; 29 : 637-646. 

[5] D. lesan. On the thermodynamics of non-simple elastic materials with two 
temperatures, J. Appl. Math. Phys. (^ZAMP) 1970; 21 : 583-591. 

[6] W.E. Warren, P. J. Chen, Wave propagation in the two temperature theory of 
thermoelasticity, Acta Mech. 1973; 16: 21-33. 

[7] W.E. Warren, Thermoelastic wave propagation from cylindrical and spherical 
cavities in the two-temperature theory, J. Appl. Phys. 1972; 43 : 3595- 3597. 

[8] D.E. Amos, On a half-space solution of a modihed heat equation, Quart. Appl. 
Math. 1969; 27 : 359-369. 

[9] S. Chakrabarti, Thermoelastic waves in non-simple media, Pure Appl. Geo- 
phys. 1973; 109 : 1682-1692. 

[10] R. Quintanilla, On existence, structural stability, convergence and spatial 
behavior in thermoelasticity with two temperatures, Acta Mech. 2004; 168 : 
61-73. 

[11] P. Puri, P.M. Jordan, On the propagation of harmonic plane waves under the 
two-temperature theory, Int. J. Eng. Sci. 2006; 44 : 1113-1126. 

[12] H.M. Youssef, Theory of two-temperature-generalized thermoelasticity, IMA 
J. Appl. Math. 2006; 71 : 383-390. 

[13] A. Magana, R. Quintanilla, Uniqueness and growth of solution in two tem¬ 
perature generalized thermoelastic theories. Math. Mech. Solids, online 2008. 

[14] H.M. Youssef, Problem of generalized thermoelastic infinite cylindrical cavity 
subjected to a ramp-type heating and loading, Arch. Appl. Mech. 2006; 75: 
553-565. 

[15] H.M. Youssef, E.A. Al-Lehaibi, State-space approach of two-temperature gen¬ 
eralized thermoelasticity of one dimensional problem, Int. J Solids Struct. 
2007; 44 : 1550-1562. 

[16] Roushan Kumar and Santwana Mukhopadhyay, Effects of relaxation time on 
plane wave propagation in twotemperature thermoelasticity, Int. J. Eng. Sci. 
2010; 48 :128-139. 

[17] S. Mukhopadhyay, R. Kumar, Thermoelastic interaction on two temperature 
generalized thermoelasticity in an infinite medium with a cylindrical cavity, 
J. Therm. Stresses, 2009; 32: 341-360. 


18 



[18] R Kumar, Anil Kumar and S. Mukhopadhyay, An investigation on thermoe¬ 
lastic interactions under two-temperature thermoelasticity with two relax¬ 
ation parameters, Math. Mech. Solids,, DOI: 10.1177/10812865145 36429. 

[19] R. Quintanilla, On Existence, Structural Stability, Convergence and Spatial 
Behavior in Thermoelasticity with Two Temperatures , Acta Mech., 2004; 
168: 61-73. 

[20] R. Picard, A structural observation for linear material laws in classical math¬ 
ematical physics. Math. Methods Appl. Sci. 2009; 32: 1768-1803.. 

[21] R. Picard. Linear thermo-elasticity in nonsmooth media. Math. Methods Appl. 
Sci., 2005; 28: 2183-2199. 

[22] S Mukhopadhyay, R Picard, S Trostorff and M Waurick, On some models 
in linear thermo-elasticity with rational material laws", Math. Mech. Solids,, 
DOI: 10.1177/1081286514556014 

[23] M. Waurick, Homogenization in fractional elasticity, SIAM .Journal on Math¬ 
ematical Analysis 2014; 46(2):1551-1576. 

[24] M. Waurick, Homogenization of a class of linear partial differential equations.. 
Asymptotic Analysis 2013; 82: 271-294. 

[25] M. Waurick, How far away is the harmonic mean from the homogenized ma¬ 
trix?, TU Dresden, 2012. 

[26] R. Picard & D. F. McGhee. Partial Differential Equations: A unified Hilbert 
Space Approach, volume 55 of De Gruyter Expositions in Mathematics. De 
Gruyter. Berlin, New York. 518 p., 2011. 

[27] A. Kalauch, R. Picard, S. Siegmund, S. Trostorff, &: M. Waurick. A Hilbert 
space perspective on ordinary differential equations with memory term. J. 
Dyn. Differ. Equations, 2014; 26: 369-399. 

[28] R. Picard, S. Trostorff, M. Waurick, & M. Wehowski. On Non-Autonomous 
Evolutionary Problems. Journal of Evolution Equations, 2013; 13:751-776. 

[29] M. Waurick, On Non-Autonomous Integro-Differential-Algebraic Evolution¬ 
ary Problems. Math. Methods Appl. Sci. , 2015; 38(4): 665-676. 

[30] S. Trostorff &: M. Wehowski. Well-posedness of Non-autonomous Evolutionary 
Inclusions. Nonlinear Analysis, 2014; 101: 47-65. 

[31] S. Trostorff. An alternative approach to well-posedness of a class of differential 
inclusions in Hilbert spaces. Nonlinear Analysis 2012; 75: 5851-5865. 

[32] C. Cattaneo. A form of heat conduction equation which eliminates the paradox 
of instantaneous propagation. Compt Rend. 1958; 247: 431-433. 

[33] P. Vernotte. Les paradoxes de la theorie continue de I’equation de la chaleur. 
Compt Rend. 1958; 246: 3154-3155. 

[34] P. Vernotte. Some possible complications in the phenomena of thermal con¬ 
duction. Compt Rend. 1961; 252: 2190-2191. 


19 



